Parameter-Free Calculation of the Solar Proton Fusion Rate in Effective Field Theory 
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Spurred by the recent complete determination of the weak currents in two-nucleon systems up to 
0(Q 3 ) in heavy-baryon chiral perturbation theory, we carry out a parameter- free calculation of the 
solar proton fusion rate in an effective field theory that combines the merits of the standard nuclear 
physics method and systematic chiral expansion. Using the tritium /3-decay rate as an input to fix 
the only unknown parameter in the effective Lagrangian, we can evaluate with drastically improved 
precision the ratio of the two-body contribution to the well established one-body contribution; the 
ratio is determined to be (0.86 ± 0.05) %. This result is essentially independent of the cutoff 
parameter for a wide range of its variation (500 MeV < A < 800 MeV), a feature that substantiates 
the consistency of the calculation. 

PACS number: 12.39.Fe 24.85.+p 26.65. +t 



The principal role of effective field theories (EFT in 
short) in nuclear physics is two-fold. One is to describe 
nuclear dynamics starting from a "first-principle" theory 
anchored on QCD. As is well known, the standard nuclear 
physics approach (SNPA) based on phenomenological po- 
tentials has been enormously successful [Q; meanwhile, 
there is growing interest in establishing the foundation 
of SNPA, identifying nuclear physics as a bona fide el- 
ement of the fundamental Standard Model. The sec- 
ond role of EFT in nuclear physics is to make precise 
model-independent predictions for nuclear observables 
with quantitative estimates of uncertainties attached to 
them. This second goal is particularly important in pro- 
viding nuclear physics input needed in astrophysics. In 
focusing on the second objective, it has been emphasized 
in a series of recent articles |2| that a very promising ap- 
proach is to combine the highly developed SNPA with an 
EFT based on chiral dynamics of QCD. This approach 
is intended to take full advantage of the extremely high 
accuracy of the wave functions achieved in SNPA while 
securing a good control of the transition operators via 
systematic chiral expansion. (See Ref. || for an an- 
other EFT calculation on pp fusion.) Such an approach- 
which is close in spirit to Weinberg's original scheme Q 
based on the chiral expansion of "irreducible terms" -has 
been found to have an amazing predictive power for the 
n + p — ► d + 7 process [||(| . 

In this paper we apply the same strategy to the solar 
proton fusion process 

p + p^ d + e + + v e . (1) 

The process (||) was previously analyzed in Ref. @ (here- 



after referred to as PKMR98) by four of the authors in 
heavy-baryon chiral perturbation theory (HBxPT). The 
calculation was carried out up to 0(Q 3 ) in chiral order, 
viz., next-to-next-to-next-to-leading order (N 3 LO). At 
N 3 LO, two-body meson-exchange currents (MEC) begin 
to contribute, and there appears one unknown parame- 
ter in the chiral Lagrangian contributing to the MEC. 
This unknown constant, called d R in J?j], represents the 
strength of a contact interaction. One intuitively expects 
that zero-ranged terms of this sort are suppressed by 
hard-core correlations in the wave functions, and hence 
it is a common practice to drop their contributions alto- 
gether. This approximation-referred to as the hard-core 
cutoff scheme (HCCS)-can indeed be justified in cases 
where the "chiral filter mechanism" associated with pion 
dominance holds || (see for details). It turns out, 
however, that the "chiral filter mechanism" does not ap- 
ply to the Gamow- Teller (GT) matrix element that gives 
the dominant contribution to the process (|). Thus there 
is no good reason to argue away the contact term con- 
tribution in this case. For practical reasons, however, 
the ci^-term contribution was simply ignored in |t| by 
invoking the HCCS and a "naturalness" argument. Let 
<$2B stand for the ratio of the contribution of the two- 
body MEC to that of the one-body current. In Q, 5^b 
for the case without the contact term is found to be 
^2B = (4.0 ± 0.5) %, where the "errors" reflect changes 
in 5zb as the hard core radius varies within the range, 
0.55 fm < re < 0.80 fm. To accommodate uncertainty 
associated with the contact term which was dropped by 
fiat, the total MEC contribution was assigned ~ 100 % 
error. With this large uncertainty attached, an EFT pre- 



diction on the pp fusion rate reported in JtJ was unable 
to corroborate or exclude the recent results of SNPA cal- 
culations ||, according to which S 2 b — 0.5 ~ 0.8 %. 

We argue here that the situation can be improved dra- 
matically. The crucial point is that exactly the same 
combination of counter terms that defines the constant 
d R enters the GT matrix elements that feature in pp fu- 
sion, tritium /3-decay, the hep process, /z-capture on a 
deuteron, and v-d scattering. This means that, if the 
value of d R can be fixed using one of these processes, 
then it is possible to make a totally parameter-free pre- 
diction for the GT matrix elements of the other pro- 
cesses. The existence of accurate experimental data for 
the tritium /3-decay rate, Tp, indeed allows us to carry 
out this program; here we are specifically interested in 
the model-independent determination of the pp fusion 
rate. The availability of extremely well tested, realis- 
tic wave functions for the ^4=3 nuclear systems enables 
us to eliminate ambiguities related to nuclear many-body 
problems. A particularly favorable aspect pertaining to 
tritium /3-decay as well as pp fusion is that both of them 
are dominated by the GT operator to a very high degree. 
In order for our result to be physically acceptable, how- 
ever, its cutoff dependence must be under control. In our 
scheme, for each value of A ~ 1/rc that defines the en- 
ergy/momentum cutoff scale of EFT, d R is determined to 
reproduce Tp; thus d R is a function of A. The premise of 
EFT is that, even if d R itself is A-dependent, physical ob- 
servables (in our case the pp-fusion rate) should be inde- 
pendent of A as required by renormalization-group invari- 
ance. We shall show below that our method indeed gives 
an essentially A-independent result, S 2 b — (0.86 ± 0.05) 
%. With this refined estimate of the two-body correc- 
tion to the well established one-body contribution, we 
are in a position to make a parameter-free prediction for 
the astrophysical S factor for pp fusion with drastically 
improved precision. 

It is worth emphasizing that the above EFT predic- 
tion for 62B is in line with the latest SNPA results ob- 
tained in Ref. |s) (and mentioned earlier). There too, 
the short range behavior of the axial MEC was con- 
strained by reproducing the GT matrix element of tri- 
tium /3-decay. The inherent model dependence of such 
a procedure within the SNPA context was shown to be 
very weak simply because at small inter-particle sepa- 
rations, where MEC contributions are largest, the pair 
wave functions in different nuclei are similar in shape 
and differ only by a scale factor ||. As a consequence, 
the ratios of GT and pp-capture matrix elements of dif- 
ferent two-body current terms are nearly the same, and 
therefore knowledge of their sum in the GT matrix ele- 
ment is sufficient to predict their sum in the pp-capture 
matrix element || . 

For the solar pp-fusion process and tritium /3-decay, 
it is sufficient to consider the limit in which the nuclear 



system receives no momentum transfer. In this limit, 
the one- and two-body axial-vector currents of the a- 
th isospin component read (in the conventions used in 
PKMR98) 
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With p = $ 1 -p 2 )/2, p l EE (pi + p{)/2, T» EE (71 X T 2 ) Q 

and similarly for a x . Eq. (|2|) includes a 1 /m 2 N correction 
term, which was ignored in PKMR98, and which is found 
to decrease the one-body matrix element by 0.13 %. In 
terms of A pp defined in Ref. |]l0| [] we have 
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6.91 



(4) 



for the central value. This should be compared with 6.93 
obtained in 0. 

The values of c's in eq.(|3|) have been determined by ir- 
N data [JlTJ] : c 3 = -3.66±0.08 and c 4 = 2.11±0.08. The 
two constants, d\ and d 2 , remain to be fixed but, thanks 
to Fermi-Dirac statistics, only one combination turns out 
to be relevant to the pp fusion processn: 



d R = di + 2d 2 



1. 2^ 1 
3 C3 + 3 C4 + 6 



(5) 



Most fortunately, as mentioned, the same combination 
enters into tritium /3-decay, the hep process, /i-capture 
on deuteron, and v-d scattering. We choose here to de- 
termine d R using the tritium /3-decay rate, T p : for which 
there are accurate experimental results. To this end, 
we calculate T p from the matrix elements of the cur- 
rent operators in Eqs.(||) and (||) evaluated for appropri- 
ate A=3 nuclear wave functions. We employ here the 
wave functions obtained in Ref. J^j using the correlated- 
hyperspherical-harmonics (CHH) method. It is obviously 



1 The subscript pp has been added here to avoid confusion 
with the cutoff parameter A. 

2 There was a sign error in the corresponding expression in 
PKMR98; as corrected here, the d-2 term should have the pos- 
itive sign. This error, however, does not affect the numerical 
results of PKMR98, since the resulting error gets absorbed 
into the definition of d R . 



important to maintain consistency between the treat- 
ments of the A=3 and A=2 systems. In our case the 
same Argonne v\s (AV18) potential Q is used (with the 
addition of the Urbana-IX three-nucleon potential fT3| ] 
for the A— 3 nuclei). Furthermore, we apply the same 
regularization method to both A=2 and A=3 systems 
to control short-range physics in a consistent manner. 
Specifically, in performing Fourier transformation to de- 
rive the r-space representation of a transition operator, 
we use the Gaussian regularization. This is equivalent to 
replacing the delta and Yukawa functions with the regu- 
lated ones, 
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The cutoff parameter A characterizes the energy- 
momentum scale of our EFT. The properly regularized 
two-body matrix elements read 



M 2 b = 
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(7) 



where u d (r) and Wd(r) are the S- and D-wave components 
of the deuteron wave function, and u pp (r) is the pp 
scattering wave (at zero relative energy). 

As mentioned, we can associate the momentum cutoff 
parameter A in Eq.(||) to the cutoff scale of our EFT. 
We let A vary within a certain reasonable range and, 
for a given value of A, we deduce the value of d R that 
reproduces Tp. With dn so determined, we calculate the 
two-body matrix element A^b- The results are given for 
three representative values of A in Table 1. The table 
indicates that, although the value of d R is sensitive to A, 
M.2B is amazingly stable against the variation of A over 
a wide range. In view of this high stability, we believe 
we are on the conservative side in adopting the estimate 
M 2 b = (0.039 - 0.044) fm. Since the leading single- 
particle term is independent of A, the total amplitude 



M = M1B+M2B, is A-independent to the same degree as 
A^2B- The A- independence of the physical quantity M, 
which is in conformity with the general tenet of EFT, 
is a crucial feature of the result in our present study. 
The relative strength of the two-body contribution as 
compared with the one-body contribution is 
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(0.86 ±0.05) %. 



(8) 



We remark that the central value of 82B here is consid- 
erably smaller than the corresponding value, &2B = 4 %, 
in PKMR98. Furthermore, the uncertainty of ±0.05 % 
in Eq.(||) is drastically smaller than the corresponding 
figure, ±4 %, in PKMR98. 

We now turn to the threshold S factor, which is a key 
input for the solar model. Adopting the most recent val- 
ues g A = 1.2670 ± 0.0035 g§ and G v = (1.14939 ± 
0.00065) x 10~ 5 GeV" 2 @, we have 



5™(0) = 3.94 




3.94 (1±0.15 %±0.10 %±e) 



(9) 



in units of 10 -25 MeV-barn. Here the first error is due 
to uncertainties in the input parameters in the one-body 
part, while the second error represents the uncertainties 
in the two-body part; e(w 0.1 %) denotes possible un- 
certainties due to higher chiral order contributions (see 
below). To make a formally rigorous assessment of e, we 
must evaluate loop corrections and higher-order counter 
terms. Although an N 4 LO calculation would not involve 
any new unknown parameters, it is a non-trivial task. 
Furthermore, loop corrections necessitate a more elabo- 
rate regularization scheme since the naive cutoff regular- 
ization used here violates chiral symmetry at loop orders. 
(This difficulty, however, is not insurmountable.) These 
formal problems set aside, it seems reasonable to assess e 
as follows. We first recall that both tritium /3-decay and 
solar pp fusion are dominated by the one-body GT ma- 
trix elements, the evaluation of which is extremely well 
controlled from the SNPA as well as EFT points of view. 
Therefore, the precision of our calculation is governed by 
the reliability of estimation of small corrections to the 



A (MeV) 


d R 


M2B (fm) 


500 


1.00 ±0.07 


0.076 - 0.035 d R ~ 0.041 ± 0.002 


600 


1.78 ±0.08 


0.097 - 0.031 d R ~ 0.042 ± 0.002 


800 


3.90 ±0.10 


0.129 - 0.022 d R ~ 0.042 ± 0.002 



TABLE I. The strength d R of the contact term and the 
two-body GT matrix element, M2B, calculated for represen- 
tative values of A. 



dominant one-body GT contribution. Now, we have seen 
that the results of the present N 3 LO calculation nicely fit 
into the picture expected from the general tenet of EFT: 
(i) the N 3 LO contributions are indeed much smaller than 
the leading order term; (ii) the physical transition am- 
plitude M. does not depend on the cutoff parameter. Al- 
though these features do not constitute a formal proof 
of the convergence of the chiral expansion used here, it 
is extremely unlikely that higher order contributions be 
so large as to completely upset the physically reasonable 
behavior observed in the N 3 LO calculation. It should 
therefore be safe to assign to e uncertainty comparable 
to the error estimate for the two-body part in Eq.(||); 
viz., e w 0.1 %. In this connection we remark that an 
axial three-body MEC contribution to the 3 H GT ma- 
trix element was calculated explicitly in SNPA |[2| and 
found to be negligible relative to the leading two-body 
mechanisms. This feature is consistent with the above 
argument since, in the context of EFT, the three-body 
MEC represents a higher-order effect subsumed in "e" in 
Eq.@. 

Apart from the notable numerical differences between 
the present work and PKMR98, it is important to point 
out that short-range physics is much better controlled 
here. In the conventional treatment of MEC, one de- 
rives the coordinate space representation of a MEC op- 
erator by applying ordinary Fourier transformation (with 
no restriction on the range of the momentum variable) to 
the amplitude obtained in momentum space; this corre- 
sponds to setting A = oo in Eq.(^). In PKMR98, where 
this familiar method is adopted, the d R term appears in 
the zero-range form, d R 5(r). Since there were no known 
methods to fix the parameter d R (except for rough es- 
timates based on the naturalness argument), PKMR98 
chose to introduce short-range repulsive correlation with 
hard-core radius rc and eliminate the d R 8{r) term by 
hand. The remaining finite-range terms were evaluated 
as functions of r C - M-2B calculated this way exhibited 
substantial r^-dependence, indicating that short-range 
physics was not well controlled; if l/rc is interpreted as 
a typical value of the momentum cutoff A in EFT, the 
results in PKMR98 imply that the theory has significant 
dependence on A. Inclusion of the d R term, with its 
strength renormalized as described here, eliminates this 
undesirable A-dependence to a satisfactory degree. 

As mentioned, the methodology employed in this work 
can be profitably used for the other related processes. 
An application to the hep process will be reported in a 
separate article [[l7| . 
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